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EFFECT OF MOLECULAR ANISOTROPY ON THE INTENSITY 
AND DEGREE OF POLARIZATION OF LIGHT 
SCATTERED FROM MODEL ATMOSPHERES 


OM P. Bahethi 

Computer Sciences Corporation 
Robert S. Fraser 

Atmospheric and Hydrospheric Applications Division, GSFC 


ABSTRACT 

Computations of the intensity, flux, degree of polarization, and the positions 
of neutral points are presented for models of the terrestrial gaseous and hazy 
atmospheres by incorporating the molecular anisotropy due to air in the 
Rayleigh scattering optical thickness and phase matrix. It is found that the 
molecular anisotropy causes significant changes in the intensity, flux and the 
degree of polarization of the scattered light. The positions of neutral points 
do not change significantly. However, if the Rayleigh scatter ing optical thick- 
ness is kept constant and the molecular anisotropy factor is included only in 
the Rayleigh phase matrix, the flux does not change, the intensity changes by 
a small amount; but the changes in the degree of polarization are still signifi- 
cant. The change in the positions of neutral points is also very small. 


v 



EFFECT OF MOLECULAR ANISOTROPY ON THE INTENSITY 
AND DEGREE OF POLARIZATION OF LIGHT 
SCATTERED FROM MODEL ATMOSPHERES 


INTRODUCTION 

In recent years there has been a renewal of interest in employing degree of 
polarization measurements (Unz 1969; Rao et al., 1973) of scattered light to 
obtain information about the sizes and indices of refraction of atmospheric 
particulates. Professor Sekera (1957) had suggested that global studies of 
atmospheric turbidity and aerosols can be made from spacecraft by making 
measurements of polarization characteristics of scattered air light. He pro- 
posed to derive aerosol parameters by comparing the measurements with a 
library of values computed from models of the atmosphere with various aerosol 
parameters. In addition, Sekera (1956) had also shown that both the degree of 
polarization and the positions of the neutral points of the skylight incident on 
the ground are related to the atmospheric turbidity. 

A large number of computational schemes are available for studying the effect 
of aerosols on the polarization characteristics of scattered sunlight (Braslau 
and Dave, 1973; Plass and Kattawar, 1972; Kattawar, Plass and Catchings, 1973; 
Hansen, 1971; Herman and Browning, 1965). Generally, it is assumed that the 
effect of molecular anisotropy on the polarization characteristics of scattered 
light is small and is partially accounted for by increasing the Rayleigh optical 
thickness according to the correct theory. Molecular anisotropy, besides 
changing the Rayleigh optical thickness, leads to a change in the scattering 
phase matrix as well. It is thus essential to incorporate the changes due to 
molecular anisotropy in the optical thickness as well as the phase matrix to 
make an accurate study of the radiative transfer problem. 

Dave (1964) has emphasized that the molecular anisotropy effect should be in- 
cluded in the radiative transfer studies of thin Rayleigh atmospheres. Schiffer 
(1972) has studied the polarization characteristics of skylight for Rayleigh at- 
mospheres by incorporating molecular anisotropy in the Rayleigh phase matrix 
only. Hansen and Hovenier (1974) have included the effect of molecular anisot- 
ropy in the determination of the atmospheric pressure at the top of the clouds 
of Venus. 

In this study we include the effect of molecular anisotropy in both the optical 
thickness and the phase matrix. Also, we analyze the consequences of including 
the anisotropy in just the phase matrix. The intensity, flux, degree of polariza- 
tion, and positions of the neutral points are computed for sunlight scattered 
from the top and the base of two model atmospheres: a Rayleigh atmosphere 
and a terrestrial atmosphere that includes aerosols. The effect of the surface 
albedo, but not of the solar zenith angle, is also studied. 
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PHASE MATRIX AND OPTICAL THICKNESS 


Chandrasekhar (1960) has derived an expression for the phase matrix for scatter- 
ing by a random aggregate of anisotropic molecules • The phase matrix R (© ) 
relates the Stokes parameters for the incident and scattered light as follows: 


I®~-L.r(0)i, (1) 

4 TT ~~ ~ 

Here I s represents the four component Stokes matrix (l|, I r s , U% and V s ) of 
the scattered light. I (I^, I r , U, V) is the Stokes matrix for the incident light. 
Chandrasekhar (1960) defines both sets of Stokes parameters I s and I. cr is 
the average scattering cross-section per air molecule and is given by 


877 ^ (m 2 - l ) 2 < 6 + Vn) . 
3 a . 4 n 2 (6 - 7 p n ) ' 


m is the index of refraction of air, N is the molecular number density, and the 
depolarization factor P n is defined as: 


Pn = 


n 


(3) 


J r 


here, 1\ and If are the intensities of light scattered perpendicularly to the 
incident beam, which is unpolarized. The scattering phase matrix R(0) is given 
by 

0 


/ 


R(0) = 


2(1 + 2 7 ) 




cos 2 0 + 7 s in 2 

y 

o 
o 

,-i 



(4) 


where 7 = p n (2 ~p n ) . The Rayleigh scattering optical thickness in the 

vertical direction above the height h is given by 


r ( a , R)(h) 


"I 


N(z) dz. 


(5) 
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From (2) and (5), we see that for a given number of molecules, the molecular 
anisotropy increases the Rayleigh scattering optical thickness by a factor 



( 6 ) 


For air p n - 0.0303 (Giicker et al., 1969). The molecular anisotropy thus in- 
creases the Rayleigh scattering optical thickness by a factor of 1.052. 

The Rubenson degree of polarization P, in the principal plane, is given by 


?=<!, -!*>/(!, + !*>■ 


(?) 


The principal plane for the plane-parallel atmospheric models that are used 
here is perpendicular to the parallel planes and contains the sun, which is 
assumed to be a point source. The degree of polarization of single scattering 
of the neutral sunlight is obtained by substituting (4) in (1) and then putting the 
resulting I® and I| in (7): 


sin 2 0(1 - y) 

1 + cos 2 ® + y( 2 + sin 2 0) 


( 8 ) 


Since y = 0.0154 for the earth T s atmosphere, it decreases the numerator and 
increases the denominator. As a result, molecular anisotropy decreases the 
degree of polarization in all directions, and the effect is largest at 0 = 90°, 
where it decreases the degree of polarization by 0.059. 

The inclination x of the plane of polarization from the £ -direction is given by 
the equation 


X ~ arctan [U/ (I^ - I f )]. 
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since U = V = 0 for the incident sunlight, we see from (1) and (4) that for single 
scattering U s =0. Hence, the plane of polarization of the scattered light is 
either parallel or perpendicular to the scattering plane and does not depend on 
the molecular anisotropy. 


MODELS AND COMPUTATIONAL DETAILS 

We have carried out computations of the polarization characteristics of light 
scattered from the following model atmospheres: 

A: Rayleigh atmosphere without molecular anisotropy. 

B: Rayleigh atmosphere with molecular anisotropy. 

C: Rayleigh atmosphere with ozone absorption, nonabsorbing haze, 
and without molecular anisotropy. 

D: Rayleigh atmosphere with ozone absorption, nonabsorbing haze, 
and molecular anisotropy. 

These atmospheres are assumed to be plane-parallel, vertically nonhomogeneous 
for models C and D, and illuminated at the top by a unidirectional monochromatic 
radiation of wavelength \ = 0.55 ^m. All orders of scattering are taken into 
account. The atmospheres are bounded on the bottom by a Lambert reflector 
with reflectivity A. 

For model A, the Rayleigh optical thickness r Cs ’ R) = 0.0924. For model B 
t (s,r) _ 0.0973. For models C and D the ozone absorption optical thickness 
= 0.0269, and the Mie scattering optical thickness of the particulates is 
T (s, m)_ 0.2979. The total scattering optical thicknesses for models C and D 
are 0.4172 and 0.4221, respectively, and differ because of molecular anisotropy 
in model D. 

In our computations, we entered the height dependence of the pressure of non- 
absorbing gases, the ozone density, and the total number density of particulates 
for thirty-five heights (0-70 km) as specified by McClatchey et al., (1971). The 
model atmosphere is divided into forty layers. In Fig. 1 the solid line curve 
shows the amount of ozone in a vertical column above a given level for models 
C and D. The total amount of ozone above sea level is 0.316 atm-cm. The 
dashed line shows the total number of particulates above a level in a vertical 
column with a cross-sectional area of one cm 2 . There are 4.91 x 10 s partic- 
ulates per cm 2 above sea level. 
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HEIGHT IN KM 


CUMULATIVE AMOUNT OF OZONE IN ATM-CM 


10 3 10' 2 10 1 10 ° 



CUMULATIVE NUMBER OF PARTICULATES PER CM 2 OF VERTICAL COLUMN 

Figure 1. An atmospheric model of the amount of ozone (solid line curve) and the 
total number of particulates in one square centimeter (dashed line) of a vertical 
column above a given level. The total ozone amount in a vertical column above 
sea level is 0.316 atm-c m and the total number of particulates is 4.91 x 10 8 cm -2 . 



For the haze in models C and D, we have used a continental distribution of the 
discontinuous type: the 'haze C" model proposed by Deirmendjian (1964). This 
model simulates the mean size distribution of particulates in continental air 

for 0.03 <, r < 0.1 fxm 


for 0.1 < r <> 5/xm 

where C is a constant; dn/dr represents the number concentration of particulates 
per cubic centimeter of air per micrometer radius. The particulate number 
density used in these computations has a strong exponential decrease up to 5 km, 
a fairly constant value until 15 km, a layer of high concentration at about 20 km, 
and for increasing heights, it falls of exponentially. A refractive index of m = 

1.5 - O.Oi is used for the aerosols. 

We have modified Dave's (1972) programs to incorporate the effect due to molecular 
anisotropy, and the results presented here employ Dave’s scheme of representing 
the phase matrix and Stokes parameters in Fourier series in the azimuth. In all 
the computations that follow, it is assumed that the flux incident on a horizontal 
surface at the top of the atmosphere is n cos 9 Q , where 0 Q is the solar zenith 
angle. 

RESULTS AND DISCUSSION 

First, we will discuss the results and comparisons of intensity, flux, degree 
of polarization, and positions of the neutral points at the top and bottom of the 
atmospheric models described above, when the effect due to molecular anisotropy 
is included in the molecular optical thickness as well as in the Rayleigh phase 
matrix. Later, we will present results for the cases when the molecular optical 
thickness is kept constant and the molecular anisotropy is included in the 
Rayleigh phase matrix only. 

a) Intensities 

Tables 1 and 2 show the intensity distribution of scattered light in the sun’s 
principal plane at the top and bottom of model atmospheres for zero ground 
reflectivity and a solar zenith angle of 60°. In the third and sixth rows of these 
tables, we have shown the relative percentage change due to molecular anisotropy 
for the Rayleigh and hazy atmospheres. For the Rayleigh atmosphere, the 
maximum change due to molecular anisotropy is 7.3 percent, whereas for the 
hazy atmosphere, the maximum change is 3.2 percent at both the top and 


masses fairly well and is given by 

fc 


dn(r) _ 
dr 


C 


— y 4 
0.1 ) 
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Table 1 

Upward intensity at the top of the atmosphere for ground reflectivity A - 0.0 and for various 
atmospheric models. B gives the percentage change from A to B, and A gives the percentage 
change from C to D. The wavelength X = 0.55 /xm and the solar zenith angle Q 0 =60°. 


Model/ 6 

0=0° 



0 = 180 

O 


0° 

20° 

o 

0 

60° 

80° 

0° 

20° 

40° 

60° 

80° 


A 

0.0218 

0.0192 


0.0428 


0.0218 

0.0293 



0.1567 

B 

0.0230 

0.0206 


0.0451 

0.1375 

0.0230 

0.0306 



0.1609 

S 

5.5 

7.3 

6.8 

5.4 

3.5 

5.5 

4.4 

3.5 

3.3 

2.7 

c 

0.0373 

0.0402 


0.1399 

0.3859 

0.0373 

0.0459 



0.1958 

D 

0.0385 

0.0415 


0.1407 

0.3822 

loggia 

0.0472 



0.1983 

A 

3.2 

3.2 

2.1 

0.6 

-1.0 


2.8 

2.2 

1.4 

1.3 


Table 2 

Downward intensity at the bottom of the atmosphere for the same conditions given in Table 1. 


0 = 180° 


0.0426 

0.0449 

5.4 

0.0729 

0.0749 

2.7 


Model / 6 

0=0° 

0° 

20° 

4^ 

O 

0 

60° 

80° 

A 

0.0217 

0.0292 

0.0422 

0.0669 

0.1545 

B 

0.0230 

0.0305 

0.0438 

0.0691 

0.1584 

8 

6.0 

4.5 

3.8 

3.3 

2.5 

c 

0.0718 

0.1426 

0.3752 

3.146 

0.9677 

D 

0.0727 

0.1429 

0.3735 

3.116 

0.9591 

A 

1.3 

0.2 

-0.4 

-1.0 

-0.9 


3.2 








































the bottom of the atmosphere. Of course, the molecular anisotropic effect on 
a hazy atmosphere depends on the haze model. The maximum changes occur 
near to where the observation angle is 90* from the sun line, which occurs at 
G = 30% <p = 0° at the top, and G = 30°, 4> - 180° at the bottom of the atmosphere. 
Tables 3 and 4 are similar to Tables 1 and 2, except for a ground reflectivity 
A = 0.50. For a Rayleigh atmosphere, we find that the maximum changes due to 
molecular anisotropy are 0.6 percent at the top and 6.3 percent at the bottom 
of the atmosphere. For the hazy atmosphere, the maximum absolute changes 
due to molecular anisotropy are 0.9 percent at the top and 2.8 percent at the 
bottom of the atmosphere. The intensity of the outward light at the top of opti- 
cally thin atmospheres lying above a bright surface is only weakly affected by 
molecular anisotropy, since most of the light comes directly from the ground, 
except near to the horizon where the path radiance is stronger than the radiance 
of the light transmitted directly from the surface. On the other hand, molecular 
anisotropy still has a significant effect on the skylight reaching the ground. 

b) Fluxes 

The flux of the diffuse light through a horizontal surface is given for the various 
models in Table 5 as a percentage of the solar flux incident on the top of the 
atmosphere. The two columns labeled top give the total upward flux at the top 
of atmosphere, and the columns labeled bottom give the downward flux of 
scattered light at the base of atmosphere. The effect of molecular anisotropy 
is to increase the fluxes because of the larger optical thickness. The effect is 
less for a hazy atmosphere than for a molecular one. The surface reflection 
decreases the effect at the top of the atmosphere but not at the bottom. 

c) Polarizations 

Figures 2 and 3 show the degree of polarization of the diffuse light as a function 
of receiver nadir and zenith angles, respectively. The usual features are 
apparent; a polarization maximum occur ing near 90° from the sun; the Babinet 
neutral point lying between the zenith and sun at the bottom, and between the 
nadir and antisolar point at the top; the Brewster neutral point appearing between 
the horizon and sun at the bottom, and between the horizon and antisolar point 
at the top; and finally, the Arago neutral point appearing at the bottom near the 
horizon, but opposite to the side containing the sun. Only two neutral points 
appear at the bottom at one time. The maximum changes in the degree of 
polarization due to molecular anisotropy are 6.9 percent for a Rayleigh atmos- 
phere (curves A and B), and 3.5 percent for a hazy atmosphere (curves C and D). 
The effect of ground reflectivity (A = 0.25) on the degree of polarization for a 
hazy atmosphere with molecular anisotropy is illustrated by curve E. 
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Table 3 

Upward intensity at the top of the atmosphere for the same conditions given in 
Table 1, except that the ground reflectivity A = 0.50. 


Model /9 



0 = 0 ° 





4 > = 180 ° 


0 ° 

— 

20 ° 

40 ° 

60 ° 

80 ° 

0 ° 

o 

O 

<M 

o 

o 

60 ° 

80 ° 

A 

0.2494 

0.2462 

0.2482 

0.2608 

0.3214 

0.2494 

0.2563 

0.2669 

0.2852 

0.3451 

B 

0.2495 

0.2464 

0.2486 

0.2616 

0.3234 

0.2495 

0.2564 

0.2672 

0.2859 

0.3469 

8 

0.0 

0.1 

0.2 

0.3 

0.6 

0.0 

0.0 

0.1 

0.2 

0.5 

C 

0.2258 

0.2269 

0.2432 

0.3024 

0.4902 

0.2258 

0.2325 

0.2495 

0.2796 

0.3001 

D 

0,2261 

0.2271 

0.2435 

0.3022 

0.4856 

0.2261 

0.2329 

0.2499 

0.2802 

0.3017 

A 

0.1 

0.1 

0.1 

- 0.1 

- 0.9 

0.1 

0.2 

0.2 

0.1 

0.5 


Table 4 

Downward intensity at the bottom of the atmosphere for the same conditions given in Table 3. 


Model / 6 



<P = 0 ° 





0 = 180 ° 


0 ° 

o 

O 

CM 

40 ° 

60 ° 

00 

o 

0 

0 ° 

o 

o 

CO 

o 

o 

60 ° 

00 

o 

o 

A 


0.0404 

0.0558 

0.0871 

0.2042 

0.0322 


tm 

0.0627 


B 

0.0340 

0.0422 

0.0579 

0.0901 

0.2100 

0.0340 



0.0659 


8 


4.5 

3.8 

3.4 

2.8 

5.6 

6.3 

5.9 

5.10 

3.4 

C 

0.0876 

0.1598 

0.3977 

3.184 

1.056 

0.0876 



0.1106 

0.2133 

D 

0.0890 


0.3965 

3.154 

1.048 



0.0747 

0.1126 


A 

1.6 


- 0.3 

- 1.0 

- 0.8 

1.6 

2.5 

2.8 

1.8 

1.1 










































Table 5 

Upward flux at the top and the diffuse downward flux at the bottom of the 
atmosphere as a percentage of the incident solar flux received at the top. 
The molecular anisotropy effect is included in the Rayleigh optical thick- 
ness as well as in the Rayleigh phase matrix. The wavelength X =0,55 
fi m, solar zenith angle d Q = 60° and surface reflectivities are A = 0.0 and 


0.5. 


Model 

Top 

Bottom 

A = 0.0 

A = 0.50 

A = 0.0 

A = 0.50 

A 

8.43 

52.31 

8.36 

12.10 

B 

8.84 

52.43 

8.76 

12.66 

8 

4.9 

0.2 

4.8 

4.6 

C 

15.95 

49.71 

34.03 

40.50 

D 

16.28 

49.87 

34.15 

40.74 

A 

2.1 

0.3 

0.4 

0.6 

8 -percentage change in going from A to B; A -percentage change 
in going from C to D. 


d) Neutral Points 


Table 6 elaborates on the positions of neutral points at the top and bottom of 
various atmospheric models for ground reflectivities of A =0.0 and 0.25. A 
comparison of the Rayleigh models A and B shows that the molecular anisotropy 
shifts the positions of the Babinet and Brewster points by about 0.1° to 0.2° 
away from the sun or antisun position for both surface reflectivities. These 
results agree with Dave ! s (1964) conclusion that molecular anisotropy has a 
weak effect on the neutral points in a Rayleigh atmosphere. 

The addition of particulates to a Rayleigh atmosphere causes the Babinet and 
Brewster neutral points at the top of atmosphere to shift by 12.2° and 9° away 
from the antisolar point, as can be seen by comparing models A and C or B 
and D. At the base of the atmosphere, the neutral points also shift away from 
the sun. The Brewster point actually disappears and is replaced by an Arago 
point. The computed result that the addition of haze to a Rayleigh atmosphere 
causes the neutral points at the bottom of the atmosphere to move away from 
the sun is contrary to observations (Sekera, 1957; Coulson, 1971). Such a 
discrepancy merits further investigation. Nevertheless, molecular anisotropy 
has a negligible effect on the neutral point positions in a hazy atmosphere. 
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DEGREE OF POLARIZATION 



>p — >pq ~ 7! NADIR ANGLE IN DEGREES 

Figure 2. The degree of polarization of light leaving the top of the atmosphere for 
various models. Wavelength k = 0.55 /zm; solar zenith 6 o - 60°; ground reflectivity 
A = 0.0. Curve E refers to a model that is the same as model D, except for A = 0.25. 



DEGREE OF POLARIZATION 



Figure 3. The degree of polarization of skylight at the bottom of the atmosphere 
for the same conditions given in Figure 2. 



Table 6 

The polar angles of the neutral points at the top and bottom of the atmosphere for ground reflectivities 
A ~ 0.0 and 0.25 for various models. The wavelength ^ = 0.55 fj-m. and the solar zenith angle &q = 60°. 
The molecular anisotropy effect is incorporated in the Rayleigh optical thickness as well as in the 
Rayleigh phase matrix. The Babinet, Brewster, and Arago points are designated by Ba, Br, and Ar, 
respectively. 



Top 

Bottom 

Model 

A = 

0.0 

A - 

0.25 


A = 0.0 



A = 0.25 



Ba 

Br 

Ba 

Br 

Ba 

Br 

Ar 

Ba 

Br 

Ar 

A 

46.7° 

75.4° 

46.7° 

75.3° 
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_ 

46.7° 

■ 


B 

46.6° 

75.6° 

46.6° 

75.5° 

46.6° 


- 

46.6° 


- 

C 

34.5° 

84.4° 

34.6° 

84.0° 

42,4° 

mm 

86.6° 

42.7° 

i i 

86.6° 

D 

34.5° 

84.4° 

34.6° 

84.0° 

42.4° 

mm 

86.6° 

42.7° 

n 

86.6° 
















RESULTS FOR PARTIAL CORRECTION DUE TO MOLECULAR ANISOTROPY 


Schiffer (1972) has obtained a solution to the radiative transfer equation by 
solving an auxiliary equation. In his models for a nonabsorbing gaseous at- 
mosphere, the molecular optical thickness is kept constant and the change due 
to molecular anisotropy is incorporated in the Rayleigh phase matrix only. 

We have repeated such calculations for all the four models presented above 
for the wavelength \ = 0.55 ^m, sun zenith & 0 = 40°, and ground reflectivities 
of 0.00, 0.25, and 0.50. The Rayleigh scattering optical thickness is r (S,R;> = 
0.0980 in all these models. 

Tables 7 and 8 show the computed results of airlight intensities at the top and 
bottom of various atmospheric models for zero ground reflectivity. In rows 
three and six of these tables, we have shown the relative percentage change 
due to molecular anisotropy for the Rayleigh atmosphere (A* and B f ) and hazy 
atmosphere (C r and D f ). The maximum change in intensity values due to 
molecular anisotropy is 1.6 percent for a Rayleigh atmosphere. Tables 9 and 10 
are computed for the same conditions of Tables 7 and 8 except for a ground 
reflectivity of A = 0.5. The maximum change in intensity due to molecular 
anisotropy is 1.0 percent for a Rayleigh atmosphere, whereas for a hazy at- 
mosphere, it is 0.4 percent. Table 11 shows that molecular anisotropy has no 
effect on the total upward and the diffuse downward flux at the top and bottom 
of the atmosphere. 

Figures 4 and 5 show the angular distribution of the degree of polarization at 
the top and bottom of models for Rayleigh (curves A f and B 1 ) and hazy (curves 
C f and D f ) atmospheres with zero ground reflectivity. Curve E T is the same 
as model D f , except that the ground reflectivity is A = 0.25. The maximum 
change due to molecular anisotropy is 6.7 percent for the Rayleigh atmosphere, 
and for the hazy atmosphere, the change is 5.1 percent at the top and 5.0 per- 
cent at the bottom of the atmosphere. 

Table 12 shows the positions of neutral points at the top and the bottom of 
model atmospheres for ground reflectivities A = 0.0 and 0.25. Molecular 
anisotropy causes a negligible change (~ 0.1°) for the Rayleigh atmosphere 
(models A 1 andB'). However, Dave (1964) found that the molecular anisotropy 
shifts the Arago point by 0.5°. For the hazy atmosphere, molecular anisotropy 
shifts the position of neutral points by about 0.4° at the top of the atmosphere, 
but there is no change at the bottom of the atmosphere. 
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Table 7 

Upward intensity at the top of the atmosphere for a ground reflectivity A = 0 and various atmospheric 
models. The wavelength k = 0.55 (im and solar zenith angle ~ 40°, The molecular anisotropy 
effect is incorporated only in the Rayleigh phase matrix. 


Model / 9 

■&- 

ii 

o 

0 

4 > = 180° 

0° 

20° 

40° 

60° 

80° 

0° 

N 

O 

0 

o 

o 

60° 

00 

o 

o 


A f 


0.0253 

0.0256 

0.0384 

0.1129 

0.0300 


0.0489 


0.1420 

B* 

E111IIS 

0.0254 

0.0260 

0.0390 

0.1133 

0.0298 


0.0481 


0.1410 

8 » 

-0.7 

0.4 

1.6 

1.6 

0.4 

-0.7 

-1.3 

-1.6 

-1.3 

-0.7 

C* 

0.0469 

0.0435 

0.0537 

0.0982 

0.2381 

0.0469 

0.0617 

0.0857 

0.1097 

0.1711 

D* 

0.0467 

0.0436 

0.0541 

0.0988 

0.2384 

0.0467 

0.0612 

0.0850 

0.1088 

0.1703 

A» 

-0.4 

0.2 

0.7 

0.6 

0.1 

-0.4 

-0.8 

-0.8 

-0.8 

-0.5 


8 percentage change in going from A 1 to B f ; A* percent change in going from C 1 to D* 


Model A 1 Rayleigh 

B 1 Rayleigh + Anisotropy 

C r Rayleigh + Ozone + Aerosol 

D* Rayleigh + Ozone + Aerosol + Anisotropy 


r C».R> = 0.0980 
r (s, R ) = 0b0980 

T Totel = 0.4225 

T"To tal „ (ftnr 























Table 8 

Downward intensity at the bottom of the atmosphere for the conditions of Table 7. 


Model / 8 

0=0° 

<t> = 180° 

0° 

20° 

40° 

60° 

o 

o 

00 

0° 

20° 

40° 

60° 

80° 

■ 

A' 


0.0377 

0.0488 

0.0684 

0.1405 

0.0300 

0.0252 


0.0382 

0.1117 

B r 


0.0372 

0.0480 

0.0675 

0.1395 

0.0298 

0,0254 


0.0388 

0.1122 

S' 

-0.7 

-1.3 

-1.6 

-1.3 

-0.7 

-0.7 

0.8 

1.6 

1.6 

0.4 

c* 

0.1489 

0.3534 

2.670 

0.6033 

0.5093 

0.1489 

0.0832 

0.0652 

0.0764 

0.1195 

D» 

0.1488 

0.3530 

2.669 

0.6025 

0.5088 

0.1488 

0.0833 

0.0652 

0.0769 

0.1197 

A* 

-0.1 

-0.1 

0.0 

-0.1 

-0.1 

-0.1 

0.1 

0.0 

0.7 

0.2 


Table 9 

Upward intensity at the top of the atmosphere for ground reflectivity A = 0.5 

and for the conditions of Table 8 . 


Model/ 8 

o 

© 

II 

0 = 180° 

0° 

20° 

40° 

60° 

80° 

0° 

20° 

40° 

60° 

D 

© 

00 


A' 

0.3880 

0.3821 

0.3785 

0.3803 

0.4063 

0.3880 

0.3947 

0.4018 

0.4106 

0.4355 

B 1 

0.3878 

0.3822 

0.3789 

0.3809 

0.4068 

0.3878 

0.3941 

0.4010 

0.4097 

0.4345 

8 f 

-0.1 

0.0 

0.1 

0.2 

0.1 

-0.1 

-0.2 

-0.2 

-0.2 

-0.2 

C 1 

0.3659 

0.3593 

0.3585 

0.3729 

0.4141 

0.3659 

0.3775 

0.3905 

0.3843 

0.3471 

D* 

0.3657 

0.3594 

0.3589 

0.3734 

0.4144 

0.3657 

0.3770 

0.3898 

0.3834 

0.3464 

A 1 

-0.1 

0.0 

0.1 

0.1 

0.1 

-0.1 

-0.1 

-0.2 

-0.2 

-0.2 














































Table 10 

Downward intensity at the bottom of the atmosphere for conditions of Table 9. 


Model/ 6 

0 = 0° 

<f> = 180° 

0° 

20° 

40° 

60° 

80° 

0° 

20° 

o 

o 

60° 

80° 

A* 

; 

0.0563 

K 

mm 

0.2225 



INMIfl 

0.0717 


B* 


0.0558 



0.2215 




0.0723 

0.1941 

S' 

-0.4 

-0.9 


-1.0 

-0.4 



0.8 

0.8 


c» 


0.3836 


0.6687 

0.6603 



0.1045 

0.1418 

0.2706 

D' 

0.1765 

0.3831 


0.6679 

0.6599 

0.1765 

1 1 

0.1046 

0.1423 

0.2708 

A' 

-0.1 

-0.1 

0.0 

-0.1 

-0.1 

-0.1 

mm 

0.1 

0.4 

0.1 





















Table 11 

Upward flux at the top and the downward diffuse flux at the bottom 
of the atmosphere as a percentage of the incident solar flux re- 
ceived at the top. The molecular anisotropy effect is incorporated 
in the Rayleigh phase matrix only. The wavelength \ = 0 .55 n m, 
the solar zenith angle 6 0 = 40°, and the surface reflectivities are 
A = 0.0 and 0.5. 


Model 

Top 

Bottom 

o 

. 

o 

II 

A = 0.50 

A = 0.0 

A = 0.50 

A* 


50.93 

5.96 

10.01 

B* 


50.93 

5.96 

10.01 

S' 

0.0 

0.0 

0.0 

0.0 

c» 

10.17 

47.40 

27.99 

35.32 

D» 

10.17 

47.40 

27.99 

35.32 

A' 

0.0 

0.0 

0.0 

0.0 


Table 12 

The polar angles of the neutral points are given at the top and bottom of the 
model atmospheres for ground reflectivities A = 0.0 and 0.25. The wave- 
length \ = 0.55 (.an. and the solar zenith angle 0 = 40°. The Rayleigh scat- 
tering optical thickness is kept constant in all these models, and the effect 
due to molecular anisotropy is included in the Rayleigh phase matrix only. 
The Babinent and Brewster points are designated by Ba and Br, respectively. 


Model 

Top 

Bottom 

A = 

0.0 

A = 

0.25 

A = 

0.0 

A = 0.25 

Ba 

Br 

Ba 

Br 

Ba 

Br 



Ba 

Br 

mm 

31.0° 

50.7° 

31.0° 

50.7° 

31.0° 

HI 

31.0° 


HR 

31.1° 

50.7° 

31.1° 

50.7° 

31.1° 


31.0° 

WjM H 


16.0° 

63.9° 


63.4° 

25.5° 

57.4° 

25.6° 

56.8° 

D 1 

15.6° 

64.3° 

15.6° 

63.9° 

25.5° 

57.4° 

25.6° 

56 .8° 
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DEGREE OF POLARIZATION 



Figure 4. The degree of polarization of light leaving the top of the atmosphere 
for various models. Wavelength A. = 0.55 (xm; solar zenith 8 o = 40°. Curves 
A', B f , C ! , D 1 and E’ are similar to those of Figure 2 except the correction due 
to molecular anisotropy is included only in the Rayleigh scattering phase matrix. 
The Rayleigh scattering optical thickness is 0.0980 in all these models. 



DEGREE OF POLARIZATION 





From these comparisons, we find that the inclusion of molecular anisotropy, in 
the Rayleigh phase matrix only, causes weaker flux and intensity changes, but 
about the same change in the degree of polarization when compared with the 
corresponding models in which the molecular anisotropy effect is included in 
both the Rayleigh optical thickness and phase matrix. 


CONCLUSION 

There are several reasons for accounting for the values of radiation parameters 
with an accuracy of a few percent or better. For example, one percent changes 
in the earth f s global albedo are believed to be important for climatology (In- 
advertent Climate Modification, SMIC). Also, our model computations show that 
if the mass of the atmospheric aerosols increases 50 percent from a nominal 
state, the intensity that would be measured by a satellite increases by 10 percent 
over a dark surface, and the maximum degree of polarization would decrease 
by a few percent over a somewhat brighter surface. Hence, the above discussion 
and comparisons make it clear that molecular anisotropy has a significant effect 
on the intensity, flux, and the degree of polarization of scattered light. The 
changes in the positions of neutral points are quite small. The changes brought 
about by incorporating the molecular anisotropy in the Rayleigh optical thickness 
as well as in the Rayleigh phase matrix are more significant compared to the 
cases when the Rayleigh optical thickness is kept constant and the change due 
to molecular anisotropy is included only in the Rayleigh phase matrix. The 
molecular anisotropy effect depends on the Rayleigh optical thickness, the 
aerosol model, and the solar zenith angle, but these relationships are not studied 
here. We conclude that to make a realistic study of radiative transfer in the 
earth's atmosphere involving scattered sunlight, the molecular anisotropy effect 
must be included in the Rayleigh optical thickness as well as in the Rayleigh 
phase matrix. 
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